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Abstract 

The relation between the bulk correlation length and the decay length 
of thermodynamic Casimir forces is investigated microscopically in two 
three-dimensional systems undergoing Bose-Einstein condensation: the 
perfect Bose gas and the imperfect mean-field Bose gas. For each of these 
systems, both lengths diverge upon approaching the corresponding con- 
densation point from the one-phase side, and are proportional to each 
other. We determine the proportionality factors and discuss their de- 
pendence on the boundary conditions. The values of the corresponding 
critical exponents for the decay length and the correlation length are the 
same, equal to 1/2 for the perfect gas, and 1 for the imperfect gas. 
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1 Introduction 



The thermal Casimir effect has been the subject of intensive theoretical [THT71 
[T9H25] and experimental research [26H29] in condensed matter physics. One 
typically considers a strongly fluctuating system contained between two parallel 
plane walls. The approach of the system's thermodynamic state towards its 
bulk critical point is accompanied by the appearance of increasingly long-range 
so-called thermodynamic Casimir forces acting between the walls. At the same 
time the vicinity of the bulk critical state induces order-parameter correlations 
of macroscopic range. It is thus natural to expect a close relation between the 
Casimir effect and the nature of these correlations in the vicinity of the critical 
state. In particular, the theory of finite-size scaling in the critical region PUS] 
points at proportionality between the bulk correlation length and the decay 
length of the Casimir forces. However, it does not specify the value of their ra- 
tio. Our purpose is to determine this ratio within the framework of microscopic 
analysis. 

We concentrate our attention on two three-dimensional systems, each un- 
dergoing Bose-Einstein condensation: the perfect Bose gas and the imperfect 
mean-field Bose gas. Our choice is dictated by the fact that a straightforward 
rigorous and microscopic analysis can be performed for each system. We ex- 
amine the relationship between the properties of the bulk correlation length, 
and the decay length characterizing the thermodynamic Casimir forces at the 
approach to the condensation point. The decay length depends on the type of 
the boundary conditions imposed at the walls while the bulk correlation length 
is independent of the boundary conditions. Near the condensation point the 
behavior of the thermodynamic Casimir force is described by the scaling func- 
tion from which one can read off the unique expression for the decay length. Its 
relation to the bulk correlation length can be established only after an indepen- 
dent calculation of the latter quantity is performed. One can th! us rephrase 
our problem as determining the bulk correlation length and the Casimir force 
decay length in two independent calculations for two systems and establishing 
the exact relation between these two quantities for each system. 

The study of the correlation length requires the knowledge of the number 
density U2(|ri — r2 1 ) of pairs of bosons situated at distance |ri — r2 1 . The density 
n 2(l r i — r 2|) tends to p 2 for |r x — r 2 | — > 00, where p is the one-particle density. 
The large distance behavior of the correlation function 

X(r) = n 2 {r) - p 2 (1) 

changes qualitatively when the system undergoes Bose-Einstein condensation. 
In the absence of condensate, the behavior of the correlation function %(r) is 
dominated by the exponential decay ~ exp(— r/£), whereas in the two-phase 
region one observes a slow non-integrable power law decay. This qualitative 
change is reflected by the fact that in the one-phase region the correlation 
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length £ governing the exponential decay tends to infinity upon approaching 
the condensation point. 

A quite similar behavior is observed for the thermodynamic Casimir force 
coupling two parallel plane walls separated by distance D. One way to evaluate 
the force is to determine, for a given temperature T and chemical potential /i, 
the excess free energy density per unit wall area 

u a (T,D,n)=u(T,D,n)-Du b (T,ii) . (2) 

Here u>(T, D, fi) is the total grand canonical free-energy density per unit wall 
area, and 0Jb{T, /a) denotes the bulk grand canonical free-energy density evalu- 
ated in the thermodynamic limit. The calculation of u> s (T, D, fj,) requires the 
thorough analysis of finite size effects under specific boundary conditions. 
By definition, the Casimir force is given by the formula 

F(T,AM) = - ^ ( ^ AM) • 0) 

Much like the correlation function, the /^-dependence of the thermodynamic 
Casimir force F(T, D, fi) in the one-phase region is dominated by the exponential 
decay ~ exp(— D/k) for D — > oo. However, its range k approaches infinity in the 
vicinity of the condensation point announcing the appearance of a new phase. 
In the presence of condensate the force F(T, D, fi) decays according to a power 
law D~ 3 . 

As already mentioned, the object of the present study is to compare the 
expressions for the correlation length £ and the range k of the thermodynamic 
Casimir force for perfect and imperfect Bose gas when the corresponding con- 
densation points are approached from the one-phase regions. 



2 Perfect Bose gas 

In the case of a perfect Bose gas the correlation function xo( r ) m the one-phase 
region (no condensate) is given by (see the excellent review [30] ) 



Xo(r) = [F(r,a,T)] ; 



(4) 



where 



A 3 F(r,«,T) = ^_exp 



-a J 



J\ 2 



(5) 



Here A = h) 'v^iro^sT is the thermal de Broglie wavelength and a = {—^/ksT. 
A particularly useful representation of series ([5]) has the form [30] 



A F(r, a,T) = - exp 
r 



X.) ^ 

E; ex P 



2 cos 



V 

-A-(s)- x 



(6) 
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with 



A ± (s) = V2^(a 2 +4 7 r 2 s 2 ) 1 / 4 



1± 



(a 2 +4tt 2 s 2 ) 1 /2_ 

It follows from Q and ([5]) that when the critical value fi Cj o = of the 
chemical potential is approached from below, and thus a — > 0, the correlation 
function xo( r ) a t large distances decays according to the asymptotic formula 

\ \ 2 , ^ / \ \ 2 / 

A \ .r , — \ / A \ IT 



A 6 X o(r)| r »A - ( ~) exp (-4^V™) = (j) cxp U|-J . (7) 

From ([7]) we find that the bulk correlation length £ of a perfect Bose gas is 
given by 

M = A^ = \{^) = 4^ A ■ (8) 

When /j, < approaches its critical value ^o,c = 0, the correlation length £o(M; T) 
diverges as Aq{— fi)~ Vo with the critical exponent i>o = 1/2, and the temperature 
independent amplitude 

A, = — %= ■ (9) 

It is our purpose here to compare the divergence of the correlation length 
£o (/Li) with that of the decay length characterizing the thermodynamic Casimir 
force ©. 

The Casimir force F(T, D, fi) for a perfect Bose gas has been evaluated by 
various methods; inter alia by Symanzik [14 via field-theoretic approach, by 
Krech and Dietrich [15] , Griineberg and Diehl u 23] , and by Zagrebnov and Mar- 
tin [9] by statistical physics methods for various boundary conditions: periodic 
(per), Dirichlet (D), and Neumann (N). The Robin boundary conditions were 
discussed by Romeo and Saharian [20], and by Diehl and Schmidt [22,24]. The 
Casimir force in the spherical model with periodic boundary conditions was 
discussed by Sachdev [TB], and Danchev [17]. In the one-phase region {fi < 0) 
the large D-dependence of the thermodynamic Casimir force is dominated by 
the exponential decay ~ exp(— D/kq), where the decay length kq for periodic 
boundary conditions is given by (see equation (21) in Ref. [9]) 




K .,per{^) = 7 \ —, T • (10) 



The decay length kq per (/i) coincides thus exactly with the correlation length 

In the case of Dirichlet and Neumann boundary conditions one obtains 

kq,d = n = Ko, P er/2- The equalities £ (m) = K o, P er(fJ-) = 2k ,d = 2k ,jv 
show clearly that the same mechanism is responsible for the appearance of the 
thermodynamic Casimir force between the walls and for building up macroscopic 
range of correlations near the condensation point. 
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3 Imperfect Bose gas 



The Hamiltonian of the imperfect Bose gas has the form 

H lmp = Ho + -—- , (11) 

where Ho is the kinetic energy operator. The term [aN 2 /2V], where a = 
J dr $(r) is a positive constant, takes into account the potential energy $(r) 
of repulsive interparticle interactions in the mean-field approximation (see Ap- 
pendix, Eq. ([37|) ). In fact, a/V > can be looked upon as the constant mean- 
field potential energy per pair of bosons. Since the mean-field theory can be 
obtained from the Kac's scaling of long-range interparticle potential [T5] it is 
worthwhile to note that the Casimir forces in the presence of long-range interac- 
tions in the spherical model were discussed by Danchev and Rudnick [19) . and 
by Danchev, Diehl, and Griineberg [21] , 

Our original derivation of the two-particle density n2 t i m p(\^i — 1*2!) for an 
imperfect Bose gas is presented in the Appendix. We show therein that this 
problem can be reduced to that of a perfect Bose gas with Hamiltonian Ho, 
and that equations (UJ) and ([5]) continue to hold provided the parameter a = 
—fi/(kBT) is replaced by 



-,rr \ v{T,n) 



(j, - ap(T, p) 
kuT 



(12) 



k B T 

The mean-field critical value of the chemical potential is positive and given by 

(13) 

where po. c denotes the critical density of a perfect Bose gas; A 3 po,e = C(3/2) = 
2.612 . This fact is a direct consequence of the implicit equation 

Pi T , M) = Po(T, v(T, (j,)) = p [T, n - ap{T, p)] (14) 

which characterizes the equilibrium state of an imperfect Bose gas ( see e.g. [3Tj). 

Consequently, the decay of the correlation function Ximp( r ) m t ne region 
P < Pimp,c is governed by the exponential law 



A 6 Ximp{f)\r^>\ 




(15) 



(compare with 0). 

In order to investigate the divergence of the correlation length £i mp upon 
approaching the condensation point we have to analyze the behavior of the 
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function v(T,p) = p — ap(T,fi) for p — > pi m p,c- The perfect gas density for 
v < is given by the series 



X 3 po(T, v) = ^ -Yj2 e*p(vq/k B T) 



9=1 



Equations (ITil) and (jTHJ) imply the following formulae 



dp . .. 
— = ( 1 + a 
op 



dp \ 1 dpo 1 



<9^ 7 9i/ a 



1- 1 + a 



dj? =il + a - 



dv I dv 2 



Using the standard notation for the Bose functions 

exp(— aq) 



9n{a) = 

9=1 



q" 



we find the following relations 

X 3 po(v) = 93/2(0) , 

5i/2(«) . 



k B TX^ - 
dv 

dv 2 



(k B T) 2 X 



: 9-1/2(0) 



The asymptotic behavior of (71/2(0;), and 5-1/2(0:) for a — > reads 

51/2(0) = \hr, 9-l/2(a) = tWt 
V o o V o 



(16) 



(17) 
(18) 

(19) 



(20) 



and these expressions permit to evaluate the second order derivative (fT5T) at the 
condensation point p = Pi mp c (v = 0) 



r 92 P r {,m dp y 3 d 2 Po 
hm 7— = hm 1 + a—— — T 

^->o dp 2 - u->o \ dv J dv 1 

X 6 k B T g-i/ 2 (a) X e k B T 
= ■ lim — = 

a 3 a-to [gi/2{ot)\ 2tt a 3 



(21) 



Thus for p pimp.c one obtains 



p - ap(T, p) 



X 6 k B T 
ilia 2 



(p pirap^c) 



(22) 



G 



It follows from the definition of the correlation length £j mj) in equation (|T5|) that 



A f k B T 

<,imp 



4 V tt[/x - ap(T, //)] 



1/2 



A ^1--^-] , (23) 



2A 2 (^, mp , c -Ai) 2C(3/2) V 

where the perfect Bose gas relation \ 3 po,c = C(3/2) has been used. 

As p, approaches fXi mPtC from below, the range of the correlation function 
£,imp diverges according to the power law 



2C(3/2) V 



(24) 



with the mean-field critical exponent Vi mp = 1. 

We now compare the above result with the analogous properties of the decay 
length Ki m p of the thermodynamic Casimir force derived by us in [32.. The 
thermodynamic Casimir force can be presented with the help of the scaling 
function T(z) [I0J|32] 

F[T ^ ] =± E [2T{z)-zT{z)\ , (25) 
where z = D/Ki mp and 

oo 

T(z) = -]T^^exp(-nz) , (26) 

n=l 

such that T(0) = -C(3)/tt, T'(0) = 0, and T(z > 1) -l/Trze" 2 . However, 
the expression for the decay length Ki mp depends on the type of boundary 
conditions imposed at the walls. For periodic boundary conditions and near the 
condensation point the decay length Ki m p,per is given by 

f 1 __ !L _y 1 (27) 



C(3/2) V M™ P ,c 

We thus find that for an imperfect Bose gas the simple relation Ki mPi p er = 
2y/n^i m p holds confirming the proportionality of these two quantities with their 
universal 2Tyir ratio near the condensation point. These two quantities diverge 
according to the same power law with the critical exponent Ui mp — 1. Simi- 
lar conclusion holds for the case of Dirichlet {nimp.u) and Neumann (Kimp.N) 
boundary conditions for which Ki mp ,D = Kimp,N = Kimp,per/2- We note that 
one could use the correlation length £j mp instead of the decay length Ki mp as 
the characteristic length rescaling the distance D in the scaling function T in 
Eq. (|25"]) . The consequence of this new choice for the scaling variable z' = D/^i mp 
is that the thermodynamic Casimir force would be now expressed in terms of the 
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new scaling function Tf,. c .(z'). Its form depends on the choice of the boundary 
conditions. 

The above results show that in the case of periodic boundary conditions the 
Casimir amplitude A defined as F(T, D, Hi mPlC )/ {ksT) — 2A/D 3 is equal 
Aj mPiPer = — C(3)/7r. The same value of the Casimir amplitude has been ob- 
tained for three-dimensional ideal Bose gas with periodic boundary conditions at 
the condensation point ^o,c = 0, sec [9, 10 . On the other hand the correspond- 
ing value for the three-dimensional spherical model (sm) with periodic boundary 
conditions is different: A smiPer = — 2£(3)/(57r), see |T7]- In the case of fluctuat- 
ing Goldstone modes of broken continuous symmetry one has A = — £(3)/(167r), 
see [TT]. The extensive discussion of Casimir amplitudes for Robin boundary 
conditions is presented in [22l[24] . 



4 Concluding comments 

We have analyzed two strongly fluctuating systems enclosed by planar walls 
for which the thermodynamic Casimir forces can be explicitly calculated: the 
perfect and the imperfect Bose gases near their condensation points. For each 
system, the decay length k characterizing the range of the exponentially de- 
caying Casimir force in the one-phase region has been evaluated and compared 
with the relevant, independently derived bulk correlation length £. The decay 
length k depends on the type of boundary conditions imposed at the walls. 
Upon approaching the condensation point of a perfect Bose gas one finds 

KQ^per — 2ko,d = 2ko.n (28) 

while the bulk correlation length £o which does not depend on boundary con- 
ditions turns out to coincide with Ko tPer - The equality £o = no,per implies that 
both quantities are divergent upon approaching the condensation point accord- 
ing to the same power law ~ (a*0,c — A*) _!/ °; where /Uo,c = 0, and uq = 1/2. 

Similar proportionality relations involving only boundary conditons depen- 
dent numerical coefficients 

^imp,per — ^^imp,D — ^^ijnp^N — ^\ / ^Cimp (^^) 

hold true in the case of the imperfect Bose gas near its condensation point at 

In order to understand the physical origin of factor 2 in Eqs (|28l29[) one has to 
go back to the structure of the energy spectrum under periodic, Dirichlet, and 
Neumann boundary conditions. The one-particle energy levels corresponding 
to the motion with momentum perpendicular to the walls under Dirichlet and 
Neumann boundary conditions for walls separated by distance D are given by 

e n' N = ^7T (f)) 2 ' ( n = !' 2 ' - for Dirichlet > and n = 0, 1,2, • • • for Neumann). 
They are equal (for n ^ 0) to the energy levels under periodic boundary con- 
ditions corresponding to distance 2D: t v ^ r = (If,) , n = 0,±1,±2, 
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It is rather straightforward to check that this fact implies the following relation 
for the Casimir excess surface energies: ujf(T,2D, /j,)\ per — 2u)f-(T, D, [J,)\d,n 
(see Eqs(43),(44), and (54) in ref.[32]). The definition of the Casimir force 

F(T, D, fj.) = — d" 3 ^gjj '^ leads then to simple equality: F per (2D) = Fd,n{D). 
In the case of exponential decay (see Eqs (l25l26p ) such an equality can hold only 
if the range of Casimir forces n under periodic boundary conditions is twice 
as big as in the case of Dirichlct and Neumann boundary conditions. We note 
that the above relation between the Casimir forces remains valid also in the 
presence of condensate (see Eqs(52, 53) in [32]). Hence, the factor 2 in Eqs ([28l 
|2"9"|) reflects the simple relation between the energy spectra. 
However, in the case of the imperfect Bose gas the value of the corresponding 
critical exponent differs form that of the perfect Bose gas, and equals Vi m p — 1 
(see ([27]) ) . This fact puts the imperfect Bose gas in a different universality 
class from the point of view of the behavior of correlations. In fact, one could 
expect here some change, as, contrary to Hq, the Hamiltonian Hi mp given in 
pip is supcrstablc 33 , 34 implying well defined thermodynamics for any value 
of the chemical potential. The question of the influence of superstability on 
thermodynamic Casimir forces has been raised in [9] , and the present study to- 
gether with previous work [32 provides a precise answer to it. Let us stress once 
more that the passage of the bulk density correlation function and the Casimir 
force from exponential decay to a power law decay at condensation point occurs 
according to the same mechanism as far as the divergence of the characteristic 
length scales is concerned. 

5 Appendix 

We derive here a simple relation between the two-particle density matrices eval- 
uated for the imperfect and perfect Bose gases. Our argument is based on the 
hierarchy equations satisfied by thermodynamic (imaginary time) Green func- 
tions [351136]. 

The one-body function is defined by 

G( 1 )(r 1 ,r 1 |r 2 ,r 2 ) = (T T [^ 1 (r 1 )^ 2 (r 2 )]) , (30) 

where < ... > denotes the grand-canonical average and 

VVxCri) = exp[n(H - M^OM^) exp[- n (H - (mN)} (31) 

V4 2 (r 2 ) - exp[T 2 (iJ - (mN)}^(t 2 ) exp[-r 2 (H - fxN)} 

are the imaginary time evolved bosonic annihilation ip( r i) an d creation i/^(r 2 ) 
operators. In Eq. (j3T))) they are ordered by the T T -operator in a chronological 
order with decreasing times from the left to the right. Their commutation 
relations read 

[V(r), ^(r')] = , [^(r), ^V)] - , [^(r), tf(i>)] - 5(r - r') . (32) 
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The system's Hamiltonian has the standard form 



H 



2m 



dr t/>t(r) Aip(r) + - / dr dr' ^(r')^^) $(r - r') ^(r)V'(r') 



(33) 

where $(r — r') denotes the interparticle potential. 

In order to derive the first hierarchy equation one applies the partial deriva- 
tive d/dri to G^)(ri, ri|r 2 , r 2 ). The result can be expressed in terms of the 
two-body Green function 

G< 2 >(n, n; r 3 , r 3 |r 2 , r 2 ; r 4 , r 4 ) = (T T [^ n (n) </v 3 (r 3 ) ^ 2 (r 2 ) <(r 4 )]> • (34) 

A straightforward calculation (see ch.7 in [35], p. 172, problem 7.2 ) yields 

^— = — AxG« + (/i - ap) G« + 5(n - r 2 ) <5(n - r 2 ) 
oti 2m 

- J dr 3 $(|r 3 - n|) G (2) (ri,Ti,r 3 ,Ti|r 2 ,r 2 ,r 3 ,Ti + ) 

-G«(r 1 ,T 1 |r 2 ,T 2 ))G«(r 3 ,r 1 |r 3 ,T 1+ )] , (35) 

where a = J dr $(?"), and ti + = limo< c ^o( T i + e )- 

In the case of an imperfect mean-field gas the potential energy per pair of 
bosons is a constant 

$(r-r') = -^, (36) 



inversely proportional to the total volume V. Inserting (131)1) into ([33f and using 
the commutation relations (1321) we find the mean field Hamiltonian 



Hi, 



2m 



dr ^(r) A-0(r) + a 



N(N — 1) 
2V 



(37) 



(38) 



For a constant potential ([36]) the last term on the right hand side of ([35 
takes the form 



where N is the particles number operator 

N = I dr ip\r)ip(r) 



G (2) (r i , ri ; r 3 , t\ | r 2 , r 2 ; r 3 , t 1+ ) 
G«(r 1 ,T 1 |r 2 ,T 2 ))G«(r 3 ,r 1 |r 3 ,T 1+ ) 



(39) 



where the integrand is a short-range integrable function. Thus, owing to the 
prefactor a/V the whole term vanishes in the thermodynamic limit V — > oo. 
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In this way we arrive at the conclusion that the one-body Green function 
Gfa£ p of an imperfect Bose gas satisfies a closed equation 

— ^ - —A lG Z + (M - ap) G« + «(n - r 2 ) 5(n - r 2 ) . (40) 



8t\ 2m 



On the other hand, in the case of a perfect Bose gas the corresponding equation 
reads ^ 

= I^AiG^ + /, G^ + Sin - r 2 ) <S(n - r 2 ) . (41) 

We thus see that the case of an imperfect gas is obtained from that of a perfect 
gas by simply replacing the chemical potential p by v = (p, — ap) in . 

Exactly the same situation persists in the case of the two-body Green func- 
tion (|34l) whose cluster decomposition 

G (2) (ri,Ti;r 3 ,T3|r 2 ,T 2 ;r4,T 4 ) = G (1) (ri, n|r 2 , r 2 )GW(r 3 ,T 3 |r 4 ,T4) (42) 

+ G (1) (ri,Ti|r 4 ,T 4 ) G (1) (r 3 ,r 3 |r 2 ,T 2 ) +G (2 ' T) (ri,ri;r 3 ,r 3 |r 2 ,r 2 ;r 4 ,r 4 ) 

is used to define the truncated function G^ 2,T \ A rather complicated structure 
of the hierarchy equation satisfied by G^ 2,T -* has been presented in Ref. [35] . 
The important property of this equation is that upon inserting the mean-field 
potential a/V into Eq.(26) of Ref. [35] and taking the thermodynamic limit one 
arrives again at a closed equation 

^^^e^-^c 1 . (43) 

The corresponding equation for the perfect Bose gas reads 



(>( "" h -A 2 G^-pG^ . (44) 



8t2 2m 

(2 T) (2) 

We thus conclude that the functions G\', and thus also , can be obtained 
by simply replacing p by (p — ap) in the corresponding perfect gas functions. 

As the two-particle number density n 2 (|ri — r 2 |) is given by equal time 
diagonal elements of the two-body Green function 



"2 



(|ri -r 2 |) = G (2) (r 1 ,T 1 ,r 2 ,T 1 |ri,T 1+ ,r 2 ,Ti+) (45) 



we also conclude that the correlation function Xim P (r) of the imperfect Bose gas 
is equal to the perfect Bose gas correlation function xo( r ) evaluated at shifted 
chemical potential v = p — ap. The analysis presented in Section III is based 
on this observation. 
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